Unwinding of a one-dimensional topological superconductor 
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We show that a topological superconductor made of four chains of superconducting spinless 
fermions characterized by four Majorana edge states can adiabatically be deformed into a triv- 
ial band insulator. To unwind this time-reversal invariant topological superconductor, interactions 
to spinful fermions are switched on along an adiabatic path. Thereby, we couple modes which be- 
long to two different representations of the time-reversal symmetry operator T with T^ = 1 and 
T^ = — 1, respectively. This observation can be understood by investigating how the relevant sym- 
metries act on the entanglement spectrum giving rise to four instead of eight different topological 
phases with Majorana edge modes. We also show that a simple level crossing of doubly and singly 
degenerate states occurs in the entanglement spectrum upon deforming the quantum state. 

PACS numbers: 



I. INTRODUCTION 

The concepts of symmetry and topology allow to clas- 
sify states of matter. A famous example is the quan- 
tum Hall effect, where the quantization of the Hall con- 
ductance is directly connected to the existence of topo- 
logical protected edge states which are robust against 
disorder and interactions. Only more recently, it was 
realizecpl'^ that similar topological states with protected 
edge modes also exist in time-reversal invariant sys- 
tems in three dimensi ons , which led to a considerable 
excitement in the fielcpl^l and intensive theoretical and 
experimentaPti^ investigation of such systems. A ma- 
jor achievement was the systematic topological classifica- 
tion for non-interacting systemJSHlsl under the assump- 
tion that only a certain set of discrete symmetries (time- 
reversal invariance and various particle-hole symmetries 
relevant, e.g., for superconductorai^ are present. 

An obvious question is to what extent this classifica- 
tion can also be used to classify interacting systems. In 
some cases such as the quantum Hall effect or the three- 
dimensional topological insulator this question can be an- 
swered as it is possible to identify measurable bulk prop- 
erties (the Hall conductance or the 9 terrrP^l^ respectively) 
which are quantized both in the presence and absence of 
interactions. For other cases the situation is less clear. 

Recently, Fidkowski and Kitaev in an insightful 
paper ^^ identified a case where a different topologi- 
cal classification was obtained for interacting and non- 
interacting systems. In one dimension, they considered 
chains of time-reversal invaria nt superconducting spin- 
less fermions (class BDJ ^^ I ^'^ ^^. In the absence of inter- 
actions, such chains are classified by an integer n € Z 
which physically describes the presence of n Majorana 
states di, i = 1, ...,n, at a given edge which all have the 
same signature under time-reversal T- Therefore the her- 
mitian quadratic coupling terms, iajaj', are forbidden as 
they are not time-reversal invariant. In contrast, inter- 
action terms of the form aia2a3a4 + H.c. are allowed. 
For n = 8 they can lead to a unique ground state at the 
edge. This fact allowed Fidkowski and Kitaev to con- 
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Figure 1: (color online) Schematic sketch of the adiabatic 
transformation from a trivial insulator (left, A = 0) to a 
superconducting state characterized by four Majorana edge 
states on each side of the open chain (right, A = 2). During 
this transformation the coupling to spinful fermions (upper 
line) is switched on. At A = 1 a pseudospin made from 4 
Majoranas forms a spin-singlet with the spin of an electron. 



struct explicitly an adiabatic path to connect the n = 8 
state to the trivial insulator. Thereby they showed that 
in the presence of interactions the classification in terms 
of Z has to be replaced by Zg. Later insightful papers 
by Turner, PoUmann and Berg''^ and also Fidkowski and 
Kitaev^ showed how these results can be understood 
by investigating how the sym metrie s of th e syst em af- 
fect the entanglement spectrurrP^IS! (gee alsJi^ESl for the 
complete classification of Id bosonic systems using pro- 
jective symmetries) obtained by calculating the eigenval- 
ues of the reduced density matrix after tracing out part 
of the system. 

The Zg topological classification of interacting spinless 
and superconducting fermions in Refs. [17118] is based on 
two symmetries, conservation of particle number parity, 
Q = (—1)^, and time-reversal symmetry, T, and three 
important equations, 7^ = 1, Q^ = 1 and [Q,?^ = 0. 
It is therefore not surprising, that if one of these condi- 
tions, T^ = 1, is relaxed, the classification scheme has to 
change. Indeed we will show that when we allow a cou- 
pling of the spinless fermions to spinful fermions (with 
T^ = ±1 for even/odd numbers of spinful fermions) it is 
possible to deform a state with 4 Majorana edge modes 
into the trivial insulator without closing the band gap, 
see Fig. [T] 



A related result has been obtained a few years ago by 
ua^l' when considering the topological stability of spin-1 
Haldane chains and of various gapped phases of spin- 
1/2 ladder models. At the edge of a gapped spin-1 Hal- 
dane chain, localized spin- 1/2 states emerge. Those are 
topologically stable, if one considers pure spin models 
with an integer total spin on each rung of the ladder. 
If, however, one remembers that in experimental system 
spin-1 degrees of freedom arise typically from the spin of 
two ferromagnetically coupled electrons, one should ask 
the question, whether the topological stability persist, 
when one considers the larger Hilbert space including 
the fermionic degrees of freedom. Here we showed by an 
explicit constructiorpi' that one can adiabatically trans- 
form the wave function of a spin-1 Haldane insulator into 
a trivial band-insulator (and Mott insulators and other 
phases of spin models) without closing the band gap by 
using this enlarged Hilbert space. As fermions even in 
an insulator can perform virtual hops from one site to 
the next, one has to include in the wave function half- 
integer instead of only integer representations of SU(2). 
This allows to unwind, for example, the spin-1 Haldane 
phase^. A general discussion of how symmetries and 
their representation lead to a topological classification of 
spin models can be found in Refs. |19|20j and [22] . 

We will start by discussing the concept of symme- 
try and adiabatic continuity in a general setup where 
a low-energy Hilbert space is embedded into some larger 
Hilbert space. Then, we will show explicitly how a trivial 
initial state can be deformed into a state with 4 Majorana 
edge modes. Finally, we will put this result in the gen- 
eral concept of classifying interacting systems in term of 
projective symmetries and the entanglement spectrum. 



II. ADIABATIC CONTINUITY AND 
CLASSIFICATION OF GROUNDSTATES 

The principle of adiabatic continuity can be used to 
classify states of matter. It can, for example, be formu- 
lated in the following way. Consider a given symmetry 
group Qs and two ground state wave functions |^o) ^nd 
l^^i) of two Hamiltonians Hq and Hf with this symmetry 
which are both defined on the Hilbert space T-ls- For sim- 
plicity we consider only systems in the thermodynamic 
limit with a unique groundstate for periodic boundary 
conditions, a finite gap in the spectrum (generalizations 
to degenerate groundstates and even gapless systems are 
possible but not considered here) and short ranged in- 
teractions. Then the two wave functions belong to the 
same topological class if a family of Hamiltonians 7^^ ex- 
ists which depend smoothly on a parameter A G [0, 1] in 
such a way that H^^^ = Hq and H^^^ = TJf . Along the 
path all symmetries in Qs have to be preserved and the 
gap of H^ has to remain finite. 

The definition given above is, however, rather restric- 
tive when one takes into account that low-energy Hilbert 
spaces are only an approximation to reality. Therefore 



it seems useful to generalize the definition of adiabatic 
continuity given above by taking into account a larger 
Hilbert space Ji in which T-lg is embedded, we consider 
H = T-Ls ®T-Le- For example the spin-1 Heisenberg model 
discussed above is originally formulated on a Hilbert 
space build from S = 1 spins. But in a real material, 
the spin 1 is formed by the Hund's rule coupling of the 
spin of two localized electrons. In this case it seems rea- 
sonable to consider for the topological classification the 
larger Hilbert space "H of the fermions rather than the 
restricted spin-1 Hilbert space 'Hs- 

Therefore we use the following generalization of the 
concept of adiabatic continuity. Two ground states. |^o) 
and l^i), of the Hamiltonians Hq and Hi defined on T-Lg 
with symmetries Qg are in the same topologically class in 
Ti. = Hs <8) He if there exist a family of gapped Hamil- 
tionians H\ which smoothly interpolate between Hq and 
Hi 



if A=0 = Hq ( 

Hx=i = Hii 



Hefi: 
Hefi 



(1) 



where H^.q is a gapped Hamiltonian in He (typically 
a simple insulator). H\ is defined on H and as above 
we require smoothness in A, short-ranged interactions, 
a unique ground state for periodic boundary conditions 
and, most importantly, all symmetries should be obeyed 
for arbitrary A. 

This last requirement, that all symmetries have to be 
obeyed, requires some care: we want to allow (as in the 
case of the spin-1 chain embedded in a fermionic Hilbert 
space discussed above) that the symmetry group Q of 
H may contain some elements (e.g. T^) which act only 
trivially (T^ = 1) on Hg- Wc require that the symme- 
tries of Q are always obeyed but Qs nevertheless differs 
from Q as the Hilbert space Hs uses, e.g., only spin-1 
representations with T^ = 1. 

At the beginning and end of the adiabatic path, A — 
0, 1, the two Hilbert spaces Hs and He are not con- 
nected. Therefore, the unique ground state wave function 
at A = (A = 1) is a product state of the ground state 
of Hq and He,o (of Hi and Heft). Note that the initial 
Hamiltonian Hefi and the corresponding wave function 
|$o), is identical to the final Hamiltonian and wave func- 
tion. Therefore the adiabatic transformation takes the 
form 



I -r . , -^ \ adiabatically , -^ , , 3- , 



(2) 



If such a transformation exist, |^o) and l^^i) are in the 
same topological sector for a given embedding in H with 
the symmetry group Q. 



III. TOPOLOGICAL PHASES OF SPINLESS 
FERMIONS 

As a warm-up we repeat the arguments leading to 
the classification of spinless superconducting fermions. 



We consider Nf chains of spinless but superconducting 
fermions described in a famous paper of Kitaev^ with 
the Hamiltionian iJ*° 



1^\ 2 \Jj,aJj + l,a''Jj,aJj + l,a 

which can be written in the form 



i?^" = 



22 {tfbj,aaj + l,a+efaj^abj,a) (3) 



.Nf 



by introducing Majorana fermions Uja, bja with a- = 
aja, &]a = bja and {aja,aja} = {bja,bja} = 2Sij using 



A. J a 



Jja ^~ Jja^ 



uja 



^yjja JjaJ^ 



(4) 



Under time reversal one therefore obtains aj^ — >■ flja 
while bja — > —bja- For tf = 0, e/ > 0, Majorana fermions 
on the same site are paired forming a trivial band insu- 
lator. In contrast, for e/ = 0, tf > Majorana fermion 
on the links are bound to each other. For an open chain 
with TV sites, j — 1, ..., N, this implies that both on the 
left side and the right side of the chain Nf Majorana 
fermions, ai_Q, and b^^ai respectively, are not bound and 
form in total 2^f zero energy states (in the absence of 
interactions) separated by a gap from all other excita- 
tions. These zero energy states remain stable for finite 
Cf < tf. Their stability in the absence of interactions is 
related to the fact that it is not possible to write down a 
quadratic term in aia which is both hermitian and time- 
reversal invariant as terms like iaiaO-ia' are odd under 
7". Therefore one can classify non-interacting systems by 
the number of Majorana edge modes, i.e. by Z. 

For an interacting systems, it is instructive to investi- 
gate for Nf = 4 the consequence of the interaction term 



Majorana modes is not topologically stable as the cou- 
pling of two pseudospins to a pseudospin-singlet allows 
for a unique ground state at the edge with no remain- 
ing degeneracies. If one assumes that the edge reflects 
the topological classification of the bulk, this observation 
'explains' the Zg classification obtained in Refs. [T6l - fT8] . 



IV. UNWINDING FOUR MAJORANA CHAINS 

We will now couple Nf = 4 chains of spinless fermions 
to spinful fermions to show that in this case one can 
unwind not only 8 but also 4 Majorana edge modes. The 
basic idea behind the following construction is that the 
coupling of a physical spin-1/2 to a pseudo spin-1/2 will 
also remove the edge mode. Furthermore, a real spin 
1/2 is not topologically protected as it can be removed 
by charge fluctuations. The combination of these two 
facts allows to construct explicitly a Hamiltonian where 
the bulk wave function interpolates smoothly between 
an ordinary insulator and a state with 4 Majorana edge 
modes of spinless fermions as in Eq. (l2| , see Fig. IT] 

We consider the following family of Hamiltonians 
parametrized by t/, e/, g, tci, ic2, U, A and J 



H 



jjes ^ 



+ X A(riia -n^b) + U X rii/s^n^pi 
J 2^ Si^ • Si^ (8) 

i.P—a,b 



^s,mt ^ _g V(aa 0^20^30^4 -I- biibi2babii + H.c.) (5) 



This Hamiltonian splits the 2^ = 4 states spanned by 
the four Majorana states an, ..., 0^4 {bn, ..., 6^4) into two 
doublets. Each of the two doublets has the properties of a 
spin-1/2. This implies that one can define for the doublet 
with the lower energy three pseudo-spin operators S^a 
i.S'ib), n^x,y,z using 



on 



l~t n ~ 



S?b = \bl<J...' (6) 



with aL 



= l(a. 



a 



lan) 



'■4 



= \{aii + 10.12) and b\^ = 

\[bi'i+ibn), b\, — \{bii+ibi2), respectively. Here one has 

to note that the operators dn^ and 6^0- can not be identi- 
fied with spinful fermions as they have the wrong trans- 
formation properties under time reversal, T. Neverthe- 
less, the pseudospin operator is a perfectly valid physical 
spin operator with the usual transformation properties 

under T, S^^ — >■ — 5," , and the conventional commuta- 
tion relation, for example, [Sf^, S^-^ = iSf^^Sij. This sim- 
ple observation immediately suggests that a state with 8 



Here H^. describes a dimerized Hubbard model with a and 
6 sites, a staggered potential A and a different hopping 
rate within the unit cell (id) and from one unit cell to 
the next (tc2)- H'^'^ couples the spin on a given site, 

s: 



03 






to the pseudo-spins of Eq. 
relevant symmetries of H are discussed in Sec. V 



6D. The 
elow. 



We consider the following path in the parameter space 
of H parametrized by < A < 2. Initially, for A = 0, the 
two subsystems are decoupled and both |*o) and |$o) 
describe trivial band insulators of the spinless and spin- 
ful fermions, respectively, with A,ey = 1 and all other 
parameters set to 0. In a first step. A, e/ are slowly 
switched off and the interaction parameters g, U and, si- 
multaneously, J are switched on. 



for < A < 1 



tf,tcl=0, 

g,U,J = A, 



e/,A = l-A (9) 
tc2 = A(l - A) 



The hopping ic2 is always finite for < A < 1 but van- 
ishes at the end of step 1, A = 1. During this first step, 
different sites i always remain disconnected as tf , ici =0 
and no entanglement builds up (see below). At the end 
of the first step, only g, U and J are finite and the a- 
Majoranas are coupled only to the a sites of the spinful 




Figure 2: (color online) Lower part of the energy spectrum of 
H as a, function of A for 1 < A < 2 (lower axes labels). The 
spectra for < A < 1 are identical (see upper axes labels) as 
E{X) — E{2 — A). The ground state is unique and the gap is 
finite. 



fermions and the 6-Majoranas to the b sites, see middle 
of Fig. [1} 

In the second step, the disconnected a and b Majoranas 
are connected again in a protocol which just reverses the 
first step up to one decisive difference: this time, a and b 
Majoranas on neighboring sites and not on the same site 
i are connected as can be seen in Fig. [T] 

forl<A<2: e/,tc2=0, i/,A = A-l (10) 

5,f/,J-2-A, t,i = (l-A)(2-A) 

Both during the first and the second step, the Hamil- 
tonian splits into blocks of just 2 spinful and 4 spinless 
fermionic sites and it is easy to diagonalize the 2^ x 2^ 
dimensional Hamiltonian numerically to show that the 
bulk gap always remains finite, see Fig. [2] 

At the end of this procedure, A = 2, all interactions 
are switched off, spinful and spinless fermions are discon- 
nected. While the wave function of the spinful fermions 
has back its original form |$o) describing a band insu- 
lator, the wave function of the spinless fermions is now 
characterized by four zero-energy Majorana states at the 
edge. Therefore, we have been able to adiabatically con- 
nect two wave functions which in the absence of the cou- 
pling to a larger Hilbert space would live in different topo- 
logical sectors. 



V. ENTANGLEMENT SPECTRUM AND 
TOPOLOGICAL STABILITY 



Instead of investigating the edge modes of a finite 
system, a very useful (and by now rather fashion- 
able) approach is to investigate the bipartite entangle- 
ment of a subsy stem and, especially, the entanglement 
spectrunp2l22ElH26]^ 'jj-^g entanglement spectrum is ob- 
tained by studying the reduced density matrix of a sub- 
system, ps = tr£;|'0)('0|, when the rest of the system. 
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Figure 3; (color online) Entanglement spectrum as a function 
of A for 1 < A < 2 for an infinite system cut into two parts. 
It evolves from an unentangled state for A < 1 (-E = In 1 = 0) 
towards a four-fold degenerate level (_B = ln4) at A = 2 re- 
flecting four Majorana edge states. The qualitatively different 
entanglement spectra are connected by a simple level cross- 
ing. Inset: The edge spectrum of an open chain also shows a 
simple level crossing of a unique ground state and a doubly 
degenerate state which merges with two other states to form 
a zero-energy four- fold degenerate state at A = 2. 



the environment is traced out. Writing ps — e^^, 
one defines the entanglement spectrum using the eigen- 
values En of H . Equivalently, one can consider the 
Schmidt decomposition of the groundstate wave function 
IV') into states of the system and environment of the form 
IV') = E„e-^"/2|n,S')|n,£;) where \n,S) and \n,E) are 
states in the system and the environment, respectively. 



A. Level crossings in the entanglement spectrum 

To calculate the entanglement spectrum, we consider 
an infinite system cut into two parts (dashed line in 
Fig. [l} . For A < 1 there is no entanglement as the Hamil- 
tonian has no matrix elements connecting the left and 
right side of the cut. Therefore there is only one non- 
vanishing eigenvalue 1 of ps and the corresponding en- 
tanglement energy is zero, _E = In 1 = 0. For A = 2, 
in contrast, the cut affects four Majorana bonds (see 
Fig. flj and the entanglement spectrum therefore consists 
of 2"^ = 4 degenerate states with entanglement energy 
_E = ln4. Fig. [3] shows how the system interpolates be- 
tween the two spectra. Starting from the 4-fold degen- 
erate, non-interacting state at A = 2, arbitrarily small 
interactions splits the four-fold degenerate state into two 
doublets when A is reduced. The degeneracy of these 
doublets is symmetry protected (see below) but when A 
is further reduced there is a simple level crossing with a 
non-degenerate state which obtains more and more prob- 
ability in the reduced density matrix (i.e., it moves closer 
to In 1=0) until it becomes the only remaining state at 



A = l. 

While the entanglement spectrum characterizing prob- 
abilities within the reduced density matrix and the spec- 
trum of edge modes are physically distinct quantities, 
they share nevertheless many similarities^^ ^^ ^^ ^21. This 
is also the case here: the inset of Fig. [3] shows the change 
of the edge spectrum as a function of A. While quan- 
titatively the entanglement and the edge spectrum are 
completely different, qualitative features, i.e., the char- 
acteristic level degeneracies and band crossings, can be 
identified. 



B. Global symmetries 

The classification of symmetry protected topological 
states relies on the precise symmetries of the system. In 
this section we will therefore discuss the relevant symme- 
tries both of the full- and the low-energy Hilbert space. 

The set of symmetry operations relevant for our dis- 
cussion is 

g = {l,T, Qc, TQc, Q, TQ, QcQ, TQcQ} (11) 

where Q = (-l)i:>= /L/.» and Qc = (-1)^-'=-'=- are 
the parity operators for the spinless and spinful fermions, 
respectively, with Q^ — Q^ — 1 and all operators com- 
mute with each other (further symmetries, e.g. inversion, 
change the classificatiorP2J) . 

While for the time-reversal invariant spinless fermions 

t T t 
one has f^^ — > fl^ , the spinor index a of spinful fermions 

t T t 
transform with the Pauli matrix a^ " — ^ "''' 



^»/3T 



«c-^; and 






T 



-IC, 



,t 
*/3f 



As a consequence the operator T^ is 1 

when acting on Hg but —1 for states in He with an odd 

number of c- fermions (as Ti = —iT). Note also that 

T ~ -J- ~ 

Sip — >■ —Si/3 and S^^ — > Sip under time reversal. 

Therefore one obtains T^ — Qc and Q can be identified 
with Z4 X Z2. 

The low-energy Hilbert space "Hs defined by the spin- 
less fermions, however, is build up from representations of 
the group Q where Qc — ^ and the remaining non-trivial 
group elements in this subspace. 



gs = {i,T,Q,TQ), 



(12) 



can be identified with Z2 x Z2 as T^ = 1 in this sub- 
space. While considering only Qg^ one obtains the Zg 
classification of interacting Majorana chains as proven in 
Ref. [16" TS^, we have shown above that four edge Majo- 
rana states arc not stable when considering the enlarged 
Hilbert space (and therefore the group Q). 



Projective symmetries and topological 
classification 



We will now show that the absence of topological pro- 
jection for n = 4 Majorana edges (and the presence of 



such a protection for n = 1,2,3) follows in a straight- 
forward way from existing clas sificat ion schemes of inter- 
acting one-dimensional systema^^ESl. In the following we 
will closely follow the arguments of Turner, PoUmann and 
BerglZl to investigate why their Zg classification scheme 
has to be modified in the presence of spinful fermions 
and why states with 0, 1, 2 and 3 Majoranas still remain 
topologically stable. To be precise, we will not try to give 
a complete classification of spinful and spinless electrons 
in one dimension. Within the spirit of the discussion in 
section [Hj we ask instead the question how interacting 
time-reversal invariant spinless fermions are classified in 
the presence of interactions if they are coupled to a trivial 
"vacuum" made out of band-insulating spinful fermions. 
We will only sketch the argumentation of Ref. [T7] and 
refer to this excellent paper for further details. 

The basic idea for the classification^^ is to investigate 
how symmetry operations affect the entanglement spec- 
trum (or more precisely, the eigenstates of the reduced 
density matrix). We define our "system" as a finite but 
very large piece cut out of an infinite system. In con- 
trast to the previous section, the entanglement there- 
fore arises now from two cuts, a left (L) and a right (R) 
cut. Turner, Pollmann and Berg^^ realized two impor- 
tant facts: (i) Each symmetry operation Si acting on p^ed 
splits naturally into two independent parts. Si ex Si^Sm, 
affecting only properties close to the left and right part. 
(ii) Classifying the possible algebras of the "splitted" 
symmetry operations SiL and Sm allows to classify the 
interacting system. To split T into two parts, one has to 
take into account that it is an antiunitary operator. This 
can, for example, be done by defining the complex conju- 
gation operator K relative to some given basis. Then T 
can be writteiP^ as T = UlUrK with unitary operators 
Ul/h acting only close to the left or right side. This also 
allows to define Tl = UlK and Tr = UrK. Defining 
also left and right parity operators, Q oc QlQr, one can 
now study the algebra of Tl/r and Ql/r- 

Consider first the case without spininl fermions, i.e. the 
situation studied in Refs. [IBHTS] where Qs defines the rel- 
evant symmetry transformations. In this case the global 
symmetries obey the three equations Q^ = 1 , T^ = 1 and 
QT = TQ. This allows for 2^ = 8 different optioniPl 
Ql/r. can either change the fermion number from even 
to odd or not, which implies that it can be either bosonic 
{QlQr = QrQl) or fermionic {QlQr = -QrQl)- Also 
Tl/r can either be bosonic or fermionic and T2 = Tr can 
either be 1 or —1. 

If we now enlarge our Hilbert space to include spinful 
fermions, then T^ = Qc is ±1 on different sectors of the 
full Hilbert space. For non-interacting systems this does 
not affect the topological classification as the different 
irreducible representations cannot couple to each other 
without breaking time-reversal symmetry: a spinful and 
a spinless fermion cannot hybridize. But in the presence 
of interaction a coupling is possible, see Eq. (|8|. As a 
consequence, the state with four Majoranas is not topo- 
logically protected any more. As shown below, that can 



directly be traced back to the fact that T^ (and there- 
fore also Tl) has always both eigenvalues ±1. States 
with 0, 1, 2 and 3 Majoranas remain, however, stable. 
They belong to different symmetry classes distinguished 
by whether Ql/r and Tl/r are fermionic or bosonic. 

The identification of the four phases can easily be 
worked out by studying how the symmetry operators act 
on chains with n = 0, 1, 2,3 Majorana edge modes (con- 
sistent with the classification table in Ref. [IZl). This 
can be done by tracking how Q and T affect the re- 
duced density matrix in the limit of uncoupled, non- 
interacting Majorana chains with tf — 1, Cf = 0. Con- 
sider, for example, the case n = 2 and a situation where 
t/ = 1,6/ = 0. In this case only the four Majorana 
operators aLi,aL2,bfa,bji2 are affected by the cut of 
the system [L/R refers to the position at the left and 
the right edge of the "system"). Defining the fermions 
a^ = \[aLi+io-L2) and 6]^ = \{i>]^i+ihi^2) allows to build 
the 4-dimensional basis |00), |01), |10), |11) describing the 
occupation of the two fermions with 0^100) — 6_r|00) = 0. 
In this space the reduced density matrix takes the form 
Pred = jl- Parity is given by Q = QlQr with bosonic 
operators Ql/r. — (—1)"^/" and n^ = ^1,0^, n^ — b\^bjj. 
Computing TQlT^^ ^ (-1)°i.°l = -Q^ shows that 
TQl = —QlT- Therefore Tl has to be fermionic. 
The explicit calculation shows that in the chosen basis 
T = cLLibR2K with a fermionic Tl =^ a,LiK where K is 

defined by K Ei ^=04 "u l«j) = E» j=o,i "ij Nj) • 

One can repeat this argument for arbitrary n: Ql is 
the product of n edge Majoranas and therefore fermionic 
for odd n and bosonic for even. Similarly, Tl is fermionic 
for n = 1,2 but bosonic for n = 0, 3. A state with 
fermionic Tl or Ql cannot be transformed smoothly into 
a state with bosonic Tl or Ql and therefore the states 
with n = 0, 1, 2, 3 belong to four different topological sec- 
tors. 

For n = 4, where both Qlir and Tl/r are bosonic, we 
have shown that an adiabatic path to n = does exist 
which uses a simple level crossing of a non-degenerate 
groundstate {T? — 1) and a Kramer's doublet (T^ = 
— 1), see Fig. [3] Such a level crossing is not possible if 
one has globally T^ = 1. As soon as one allows that Tl 
has simultaneously the eigenvalues +1 and —1, it is not 
possible to fulfill the condition that T^ = 1 for all its 
eigenstates. 

We have therefore shown, that states with 0,1,2 and 3 
Majorana modes at a given edge are topological stable in 
our system while one can deform a state with n Majorana 
modes at the edge in a state with n mod 4 Majorana 
modes. Note, however, that we have only focused on the 
low-energy Hilbert space of our subsystem while we have 
not investigated the possibility of other topological states 
involving the spinful fermions. 

VI. CONCLUSIONS 

With this paper we have added one more example to 
the (short) list of cases where interactions enforce to mod- 



ify the classification of non-interacting topological insu- 
lators and superconductors. Our simple examples can be 
used to demonstrate several points: (i) The topological 
classification of a low-energy Hilbert space can be insuf- 
ficient in the sense that in a larger Hilbert space it is 
possible to unwind topological states of the low-energy 
space. This can happen if the larger Hilbert has different 
symmetries in the following sense: the irreducible rep- 
resentations forming the low-energy Hilbert space lead 
to a different symmetry group. In our specific example, 
the key was that one of the symmetries of Ti. (T^) had 
a trivial representation (T'^—1) in "Hs and the groups Q 
in Eq. (TTll) and Qs in Eq. ( 12 ) differ. Another example 



for this are spin models embedded in fermionic Hilbert 
spaces^il. (ii) When one uses entanglement to character- 
ize a system, it is not sufficient to concentrate on the 
lowest states in the spectrum. Levels with other degen- 
eracies and symmetry properties may cross with the low 
energy levels and completely rebuild the entanglement 
properties without closing any bulk gap. To exclude this, 
one has to classify all those states in the entanglement 
spectrum which survive in the thermodynamic limit (i.e. 
for the reduced density matrix of an infinitely large sub- 
system in an infinitely large total system), (iii) Finally, 
our example demonstrates the validity power of the exist- 
ing classification schemes of interacting one-dimensional 
system^ 



IY'18; 



An important ingredient of our analysis was the pres- 
ence of spinless fermions. Building a system well de- 
scribed by spinless fermions interacting with spinful 
fermions is not simple but at least in principle possi- 
ble. Spinless fermions can, for example, be realized with 
ultracold alkaline atoms with a polarized nuclear spin. 
While polarizing the nuclear spin nominally breaks T- 
invariance, due to the closed electronic shell of alkaline 
atoms this has practically no consequences for the inter- 
action with other particles. 

An important challenge for the future is to find a com- 
plete classification scheme for interacting spinful elec- 
trons in higher dimensions. It will also be useful to search 
for higher-dimensional model systems which can be used 
to study in concrete examples how interactions modify 
the topological classification, see for example Ref. |27) . 
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